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ABSTRACT 



Context. Strongly star-forming galaxies of subsolar metallicities are typical of the high-redshift universe. Here we therefore provide 
accurate data for two low-z analogs, the well-known low-metallicity emission-line galaxies Haro 1 1 and ESO 338-IG 004. 
Aims. Our main goal is to derive their spectroscopic properties and to examine whether a previously reported near-infrared (NIR) 
excess in Haro 1 1 can be confirmed. 

Methods. On the basis of Very Large Telescope/X-shooter spectroscopic observations in the wavelength range ~AA3000 - 24000, we 
use standard direct methods to derive physical conditions and element abundances. Furthermore, we use X-shooter data together with 
Spitzer observations in the mid-infrared range to attempt to find hidden star formation. 

Results. We derive interstellar oxygen abundances of 12 + log O/H = 8.33 + 0.01, 8.10 ± 0.04, and 7.89 ± 0.01 in the two H II 
regions B and C of Haro 1 1 and in ESO 338-IG 004, respectively. The observed fluxes of the hydrogen lines correspond to the 
theoretical recombination values after correction for extinction with a single value of the extinction coefficient C(H/3) across the 
entire wavelength range from the near-ultraviolet to the NIR and mid-infrared for each of the studied H II regions. Thus, we confirm 
our previous findings obtained for several low-metallicity emission-line galaxies (Mrk 59, II Zw 40, Mrk 71, Mrk 996, SBS 0335- 
052E, PHL 293B, and GRB HG 031203) that the extinction coefficient C(H/3) is not higher in the NIR than in the optical range and 
therefore that there are no emission-line regions contributing to the line emission in the NIR range, which are hidden in the optical 
range. The agreement between the extinction-corrected and CLOUDY-predicted fluxes implies that a H II region model including 
only stellar photoionisation is able to account for the observed fluxes, in both the optical and NIR ranges. No additional excitation 
mechanism such as shocks from stellar winds and supernova remnants is needed. All observed spectral energy distributions (SEDs) 
can be reproduced quite well across the whole wavelength range by model SEDs except for Haro 11B, where there is a continuum 
flux excess at wavelengths >1.6/jm. It is possible that one or more red supergiant stars are responsible for the NIR flux excess in Haro 
1 IB. We find evidence of a luminous blue variable (LBV) star in Haro 1 1C. 

Key words, galaxies: fundamental parameters - galaxies: starburst - galaxies: ISM - galaxies: abundances - stars: activity 



1. Introduction 

Haro 11 (=ESO 0350-IG 038) and ESO 338-IG 004 (=Tololo 
1924-416) are prominent blue compact galax ies (BCG) with Mb 
= - 2 0.0 mag and -18.9 mag, respectively (Bergva ll & Ostlinl 
2002). Both exhibit perturbed morphologies, high star-formation 
rates (SFRs), and multiple-componen t Ha velocities f ields re- 
flecting signatures of merger events (IQstlin et all l2001). 

Haro 11 i s a massive and luminous BCG with a S FR ~18- 
20 M Q yr _1 dBergvall & Ostlinl |2002|) which possesses multi- 
ple H I I regions, ca l led A , B ; and C fo ll owing the nomencla- 
ture of iKunth etaTJ d2003l) . IQstlin et alj d2001h estimated that 
the stellar mass of the galaxy is 10 10 M Q from the Ha veloc- 
ity field, while Bergvall et al. (2000) derived a total gas mass of 
2xl0 9 M Q . The high present SFR o f 22+3 M yr _I in Haro 11 
was derived by lAdamo et al. (l2010t). The galaxy is a Lya and 
Ly-continuum emitter (lHaves et all 120071) and also a luminous 



IRAS source (Bergvall et al., 2000) with an IR luminosity of 
1.9xlO n Lq, which is a characteristic property of luminous IR 
galaxies (LIRGs). High star-formation activity in Haro 1 1 and 
its global properties are similar to those in high-redshift Lyman 
break galaxies (LBGs) at z ~ 3 dPettini et all 1200 lb and those in 
lower-redshift lu minous compact emis sion-line galaxies (LCGs) 
from SPSS DR7 (Izoto v et all 1201 la) and in green pea galaxies 
dCardamone et alll2009l) ~ 
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* Based on observations collected at the European Southern 
Observatory, Chile, ESO programme 60A-9433(A). 
** Figures 1, 2 and 3 and Table 1 are available in electronic form at 
http://www.aanda.org 



ESO 338-IG 004 contains two H II regions in its most ac- 
tive part, which is referred to its "centr e" and "H II region" 
dBergvalll [1981 iBergvall & Ostlinl |2002|) . The central starburst 
("centre") of the galaxy is resolved by Hu bble Space Telescope 
(HST) images into compact star clusters dMeurer et all 1995. 
Ostlin et al., 1998) and young globular clusters (Ostlin et al., 
1998). Bergvall (1985) estimated the dyn amical mass of the 
galaxy, M tot ~ lxl0 9 M o . Ostl in et all d2007l) derived the dynami- 
cal mass M = 1.3xl0 7 M o and an age of ~6 Myr for the brightest 
cluster #23 in ESO 338-IG 004 from an analysis of the absorp- 
tion components of H<5 and H8 - HI 1 Balmer lines. 

IBergvall & Ostl in (2002) obtained a very red colour V - 
K=4.2±0.8 ma g for a Haro 1 1 unde rlying LSB component. On 
the other hand, Mic heva et all (1201 Oh derived a far less red V - K 
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colour of the Haro 11 LSB component of 2. 3 ±0.2 mag from new 
deep V and K photometry. On the basis of high-resolution HST 
imaging of H aro 1 1 in eight wavele ngth bands from the ultravio- 
let (UV) to K. lAdamo et all (1201 Ol) found about 100 star clusters 
with flux excesses at wavelengths > 8000 with respect to the 
synthetic evolutionary models. This is almost half of the whole 
cluster sample in Haro 1 1 . Moreover, both of the giant H II re- 
gions of Haro 11B and C are very massive (~lxl0 7 M ) and 
have a red excess. 

Based on ESO 3.6m telescope observations, Bergvall ( 1985) 
obtained an oxygen abundance of 12 + log O/H = 8.08 for the 
central 3"x4" region of ESO 338-IG 004 by applying the di- 
rect r e -method with the use of the emission line [O Ill]/14363. 
An oxygen abundance 12 + log O/H = 7.92, which was d erived 
usin g the 7V-method, was obtained by iMasegosa et al.1 (1 19941) 
and Raimann et al. (2000) from the same spectra acquired by 
Terle vichetalJ (119911) . 

Finally, iBergvall & Ostlinl d2002l) using ESO 1.5-m and 3.6- 
m spectra (obtained in 1983, 1984 and 1986) derived 12 + log 
O/H = 7.9 and 8.0 for Haro 11 and ESO 338-IG 004, respec- 
tively, again using the r e -method. 

We note, however, that in the case of Haro 11, such a funda- 
mental parameter as the metallicity remains uncertain despite of 
many previous comprehensive studies. This is because the im- 
portant [O III]/14363 line is blended with the [Fe ll]/14359 line 
(this paper) and therefore an accurate measurement of the former 
line's flux was impossible with previous low-spectral-resolution 
observations. 

We present here high-quality archival VLT/X-shooter spec- 
troscopic observations of the two bright knots B and C in 
Haro 1 1 and of the brightest part of ESO 338-IG 004 over a 
wide wavelength range ~/l3000 - 24000. These new medium- 
spectral-resolution observations allow us to more accurately de- 
rive reddening, the physical conditions, and the element abun- 
dances in the H II regions. Moreover, since cool low-mass stars 
are a main contributor to the stellar mass and they emit mainly 
in the near-infrared (NIR) range, observations of a very broad 
wavelength range permit us to derive more reliable stellar masses 
for the galaxies and their specific SFRs. The observations from 
the UV to NIR wavelength ranges also allow us to study whether 
any hidden star formation is present in these galaxies. We note, 
however, that in low-metallicity compact galaxies with high 
S FRs the strongly ionised gaseous emission contributes to the 
total spectral energy distribution (SED) and therefore should be 
properly taken into account for the stellar mass determination. 

Thus, the main goals of this paper are to obtain more accu- 
rate element abundances, to test spectro scopically the flux excess 
in the NIR wavelength range found by lAdamo et al.l d2010b for 
Haro 11B and C, to search for star formation seen in the NIR 
range but hidden in the optical range, and to estimate the ex- 
citation mechanisms of brigh t emission line s using additionally 
the Spitzer observations of IWu et all (12008) in the mid-infrared 
(MIR) range. 

In Sect. |2 we describe the observations and data reduc- 
tion. We present the main properties of the galaxies in Sect. [3] 
More specifically, the element abundances are presented in sub- 
sect. 13.11 the hidden star formation is discussed in subsect. 13.21 
the H2 emission is considered in subsect. 13.31 the CLOUDY 
modelling is discussed in subsect. 13.41 the possible presence of 
LBV stars in Haro 11C and Haro 11B in subsect. 13.51 and the 
SED fitting is described in subsect. 13.61 Finally, in Sect. [4] we 
summarize our main results. 



2. Observations 

New spectra of the blue compact emission-line galaxies Haro 
1 1 and ESO 338-IG 004 were obtained with the VLT/X-shooter 
on 2009 August 11 and 12, respectively [ESO program 60. A- 
9433(A)]. 

The observations were performed in the wavelength range 
~/l3000 - 24000 using the three-arm echelle X-shooter spectro- 
graph mounted at the UT2 Cassegrain focus. All observations 
were obtained at low airmass 1.016, 1.012, and 1.430 for Haro 
1 IB, Haro 1 1C, and ESO 338-IG 004, respectively, thus the ef- 
fect of the atmospheric dispersion was low. The seeings were 
0772 - 0775, 0770 - 1732, and a/61 - 0766 for Haro 11B, C, 
and ESO 338-IG 004 observations, respectively. For each of 
the UVB, VIS, and NIR arms, the total exposure times were 
800s, 680s, and 600s for Haro 11B, Haro 11C, and ESO 338- 
IG 004, respectively. They were broken into two equal subex- 
posures. The nodding along the slit was performed according to 
the scheme AB with the object positions A or B differing by 4" 
along the slit for each H II region. In the UVB arm with wave- 
length range ~/13000 - 5600, a slit of l"x 11" was used. In the 
VIS and NIR arms with wavelength ranges ~/15500 - 10200 and 
~il0200 - 24000, respectively, slits of 079 x 11" were used. 
The binning factors in the UVB and VIS arms along the spa- 
tial and dispersion axes were 1 and 2, respectively, while in the 
NIR arm they were 1 and 1. Resolving powers A/AA of 5100, 
8800, and 5100 characterise the UVB, VIS, and NIR arms, re- 
spectively. Spectra of thorium-argon (Th-Ar) comparison arcs 
were used to perform the wavelength calibration of the UVB 
and VIS arm observations. For the wavelength calibration of the 
NIR spectrum, we used night sky emission lines. 

The two-dimensional UVB and VIS spectra were bias sub- 
tracted and flat-field corrected using IRABJ. The NIR dark cur- 
rent frames were subtracted from the NIR flat frames. The two- 
dimensional NIR spectra of objects were then divided by the 
flat frames to correct for the pixel sensitivity variations. Cosmic 
ray hits of all UVB, VIS, and NIR spectra were removed using 
the routine CRMEDIAN. The remaining hits were later removed 
manually after background subtraction. 

For each of the UVB, VIS, and NIR arms, the spectrum at 
the position B was subtracted from the spectrum at the posi- 
tion A. This resulted in a frame with subtracted background. 
We used the IRAF software routines IDENTIFY, REIDENTIFY, 
FITCOORD, and TRANSFORM to perform wavelength calibra- 
tion and correct for distortion and tilt for each frame. The one- 
dimensional wavelength-calibrated spectra were then extracted 
from the two-dimensional frames using the APALL routine. We 
adopted extraction apertures of l"x 4", 079 x 4", and 079 x 4" 
for the UVB, VIS, and NIR spectra, respectively. Before extrac- 
tion, the spectra at the positions A and B in the two-dimensional 
background-subtracted frames were carefully aligned with the 
routine ROTATE and co-added. The spectrophotometric stan- 
dard star Feige 110 was used to perform the flux calibration. 
Its spectra in all arms were obtained with a wide slit of 5 "x 11". 

The resulting flux-calibrated and redshift-corrected UVB, 
VIS, and NIR spectra of Haro 1 IB, C, and ESO 338-IG 004 are 
shown in Figs.Q~|[2l and [3](available only in the online edition). 



1 IRAF is the Image Reduction and Analysis Facility distributed by 
the National Optical Astronomy Observatory, which is operated by the 
Association of Universities for Research in Astronomy (AURA) under 
cooperative agreement with the National Science Foundation (NSF). 
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Fig. 4. Close-up part of the Haro 1 IB optical spectrum shown in Fig. Q] 



3. Results 

3.1. Element abundances 

We derived element abundances from emission-line fluxes us- 
ing a classical semi-empirical method. The fluxes in all spec- 
tra were measured using the IRAF SPLOT routine. The line 
flux errors include statistical errors derived with SPLOT from 
non-fiux-calibrated spectra, in addition to errors introduced by 
the absolute flux calibration, which we set to 1% of the line 
fluxes, according to the uncertainties in the a bsolute fluxes of rel- 
atively bright standard stars (jo ke, 1990; Co lina& Bohlinlll994t 
lBohlinL[l996l : llzotov & T huan. 2004). These errors were propa- 
gated into the calculation of the electron temperatures, the elec- 
tron number densities, and the ionic and total element abun- 
dances following the prescription of G useva et al.l (1201 lbh . 

The extinction coefficient C(H/3) and equivalent widths of 
the hydrogen absorption lines EW(abs) are calculated by si- 
multaneously minimising the deviations in the corrected fluxes 
I{A)/I{Hf3) of all hydrogen Balmer lines from their theoretical 
recombination values. The fluxes were corrected for both ex- 
tinction, using the reddening curve of Cardelli et al. ( 1989), and 
underlying hydrogen stellar absorption (llzotov et al.l 09941) . The 



derived C(H/3) was applied to correct all emission-line fluxes in 
the entire wavelength range AA3QQQ - 24000. The extinction- 
corrected relative fluxes I(A)/I(H/3)xlQ0 of the lines, the extinc- 
tion coefficient C(H/?), the equivalent width of the H/3 emission 
line EW(Hy6), the H/3 observed flux F(HyS), and the equivalent 
width of the underlying hydrogen absorption lines EW(abs) are 
given in Table Q] (available only in the online edition). We note 
that the fluxes of the Balmer hydrogen emission lines corrected 
for extinction and underlying hydrogen absorption (Cols., 3, 5, 
and 7 in Table [D are close to th e theoretical recombination val- 
ues of lHummer& Storevl(ll987l) (column 8 of the Table 1), sug- 
gesting that the extinction coefficient C(HyS) is reliably derived. 

The physical conditions, and the ionic and total heavy ele- 
ment a bundances in the IS M of the galaxies were derived fol- 
lowing Hzotovet al. (2006) (Tabled. In particular, applying the 
direct r e -method, we derived the electron temperature r e (0 III) 
from the [O III] ^4363/(/14959 + ,15007) emission-line ratio. The 
electron temperatures r e (0 II) and T e (S III) were derived from 
the empirical r elations based on the photoionisation models of 
H II regions (I zotov et all 120061) . The electron number densities 
N e (S II) were obtained from the [S ll]/16717//16731 emission- 
line ratios. The electron temperatures r e (0 III), T e (Q II), and 
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r e (S III), the electron number densities iV e (S II), the ionisation 
correction factors ICFs, and the ionic and total O, N, Ne, S, Ar, 
Fe, and CI abundances derived from the forbidden emission lines 
are given in Table [2] 

We derive 12+logO/H = 8.33+0.01 and 8.10+0.04 for 
Haro 11B and C, respectively, which is appreciably higher 
than the oxygen abund ance 12+logO/H = 7.9 obtained by 
Berg vall & Ostlinl (120021) in an aperture 4"x4" for the bright- 
est region of Haro 11. We also obtain a higher extinction 
C(HB)=QJ4 for Haro 11B than C(HyS)=0.27 for Haro 11C. It 
can in deed be seen in the HST images obtained by Adam o et alJ 
(2010) tha t Haro 1 IB is m ore obscured by dust clouds than Haro 
11C. lAdamo et alJ (1201 Oh also derived different extinctions of 
E(B-V)= 0.38 and 0.06 (or C(HB) = 0.56 and 0.09) for Ha ro 1 IB 
and Haro 1 1C, respectively. For comparison, Ha ves et all (12007) 
derived E(B- V)= 0.42 and 0.48 (or C(HjS)=0.62 and 0.70) and 
IVaderetal.l(ll993l) derived E(B-V)=QA1 and 0.39 for Haro 1 IB 
and C, respectively. 

The differences in oxygen abundance of Haro 1 1 are likely 
not solely due to differences in the extinction but may also 
be caused by blending of the [O Ill]/14363 emission line with 
the [Fe ll]/14359 one in previous low-spectral-resolution obser- 
vations. The higher-spectral-resolution X-shooter observations 
have allowed to resolve this blend. In TableQ~]we can clearly see, 
that the fluxes of the [Fe Il]/I4359 emission line in Haro 1 IB and 
Haro 11C are, respectively, 27% and 151% of the [O Ill]i4363 
emission-line fluxes. 

The oxygen abundance 12+logO/H = 7.89+0.01 of ESO 
338-IG 004 is slightly l ower than the previous determinations 
of 12 + log O/H = 7.92 dMasegosa et all 1 1994 IRaimann et all 



2000), 8.08 (Bergvall, 1985), 8.0 dBergvall & Ostlinl 120021) and 
7.99 dGuseva et all 1201 lal) . where the direct r e -method was 
also used. We note that our determination of C(HB) = 0.31 (see 
Table [B is lar ger th a n the C(HB) in t h e rang e 0.0 - 0.05 derived 



by Bergval i 1985b . IRaimann et~al 
(2011a)", and E(B 



_ (200 0J), and iGuseya et all 
V)<0.1 of Bergva ll & Ostlinl (|2002|) . From 



the fits of model SEDs to the SEDs derived from HST images 
with different filters, Ostlin et al.1 d2009l) derived a low value 
of C(H/J)~0.1 for ESO 338-IG 004. These differences may be 
caused by the oxygen abundances having been derived by differ- 
ent authors for different regions in ESO 338-IG 004. To check 
this possibility, we extracted separate spectra in addition to the 
integrated spectrum, for the "centre" and "H II region". For the 
centre, we derived 12+log O/H = 7.82+0.01, which is lower than 
that for the integrated spectrum, but for the H II region we obtain 
a higher value of 8.00+0.01, which is more consistent with some 
of the previous observations. 

The abundance ratios of the different species log N/O, log 
Ne/O, log S/O, log Ar/O, log Fe/O, and log Cl/O (Table 
obtained for Haro 11B, C, and ESO 338-IG 004 are in the 
range of the ratios obtained for different sampl es of blue 
compact g alaxies, such as the H eBCD sample of llzotov et alJ 
(2004) and Izot ov & Thuanl (12004) used to determine the primor- 
dial He abundance (consisting of more than 100 emission-line 
galaxies), the sample of ~300 emission-line galaxies from the 
SDSS DR3 of Izot ov et alJ ( 120061) . and the VLT s ample of 121 
low-m etallicity emission-line galaxies studied by iGuseva et alJ 
(1201 lal) . The high values of log N/O in Haro 1 IB (-0.92, Table 
O and in Har o 1 1C (-0.79) agree with the value of log N/O =- 
0.7 derived by [Bergvall & Ostlinl d2002l) for Haro 1 1 . Our values 
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Table 2. Physical conditions and element abundances 



Table 3. Stellar masses and ages of young stellar populations 



Galaxy 



Haro 1 IB 



Haro 11C 



ESQ 338-IG 004 



Property 



Value 



r e (0 III), K 
r c (0 II), K 

r e (s in), k 

iV e (S II), cm-3 

O + /H + ,(xl0 5 ) 
O 2+ /H + ,(xl0 5 ) 
3+ /H + , (xlO 5 ) 
O/H, (xlO 5 ) 
12+log O/H 

N + /H + , (xlO 6 ) 
7CF(N) 11 
N/H, (xlO 6 ) 
log N/O 

Ne 2+ /H + , (xlO 5 ) 
7CF(Ne) a 
Ne/H, (xlO 5 ) 
log Ne/O 

S + /H + , (xlO 6 ) 
S 2+ /H + , (xlO 6 ) 
7CF(S) a 
S/H, (xlO 6 ) 
log S/O 

Ar 2+ /H + ,(xl0 7 ) 
Ar 3+ /H + ,(xl0 7 ) 
7CF(Ar) a 
Ar/H, (xlO 7 ) 
log Ar/O 

[Fe III] ,14658: 
Fe 2+ /H + ,(xl0 6 ) 
ICF(Fef 
Fe/H, (xlO 6 ) 
log Fe/O 

[Fe III] ,14986: 
Fe 2+ /H + ,(xl0 6 ) 
ICF(Fef 
Fe/H, (xlO 6 ) 
log Fe/O 

Cl 2+ /H + ,(xl0 7 ) 
lCF(C\f 
Cl/H, (xlO 7 ) 
log Cl/O 



10144±70 

10099±77 

10110±58 

10+10 

8.38±0.28 
12.78±0.33 
0.22±0.02 
21.38±0.43 
8.33±0.01 

8.87±0.29 

2.93 
25.99±1.06 
-0.92±0.02 

3.10+0.09 

1.49 

4.62+0.17 

-0.66+0.02 

0.57+0.01 
2.39+0.08 

1.07 

3.17+0.10 

-1.83+0.02 

6.23+0.14 
0.58+0.05 

1.07 

6.69+0.16 

-2.50+0.01 



0.99+0.03 

3.73 

3.69+0.16 

-1.76+0.02 



0.35+0.02 

3.73 

1.30+0.07 

-2.21+0.02 

0.35+0.02 

1.13 

0.39+0.02 

-3.73+0.03 



11362+449 

11552+500 

11209+373 

10+10 



5.21= 
7.48= 



0.76 
0.91 



12.69+1.19 
8.10+0.04 

7.48+1.13 

2.72 
20.37+3.94 
-0.79+0.09 

2.09+0.28 

1.42 

2.97+0.49 

-0.63+0.08 

0.68+0.06 

1.48+0.31 

1.06 

2.29+0.35 

-1.74+0.08 

5.59+0.44 

1.08 

6.02+0.52 

-2.32+0.06 



1.32+0.17 

3.48 

4.61+0.88 

-1.44+0.09 



0.51+0.09 

3.48 

1.78+0.42 

-1.85+0.11 

0.71+0.12 

1.17 

0.84+0.23 

-3.18+0.13 



14771+137 

13786+160 

14397+114 

142+36 

1.66+0.06 
5.87+0.16 
0.19+0.01 
7.73+0.17 
7.89+0.01 

0.62+0.03 

4.54 

2.82+0.16 

-1.44+0.03 

1.34+0.04 

1.13 

1.52+0.05 

-0.71+0.02 

0.22+0.01 
0.77+0.04 

1.23 

1.22+0.05 

-1.80+0.02 

2.42+0.06 
0.89+0.07 

1.11 

2.69+0.11 

-2.46+0.02 



0.24+0.01 

6.17 

1.49+0.10 

-1.71+0.03 



0.29+0.01 

6.17 

1.77+0.11 

-1.64+0.03 



ID 



Mass(M G ) Age(Myr) C(H/J) 



Haro 1 IB 
Haro 11C 
ESO 338-IG 004 

Haro 11B 
Haro 11C 

ESO 338-IG 004 



1.26xlO B 
2.05xl0 8 
0.56x10 s 



this paper" 



2.9 

4.2 
1.9 



0.74 
0.27 
0.37 



Adamoetal. (2010)° 



8.35x10° 
1.36xl0 7 



3.5 

9.5 



0.56 
0.09 



Ostlin et al. (2007) c 



1.3x10' 



6.3 



Ionisation correction factor. 



a VLT/X-shooter data. 

b HST photometric data. 

C VLT/UVES data for the brightest cluster #23. 



enables us to directly compare the extinctions derived for the 
optical and NIR spectra. 

The observed F(A) and extinction-corrected 1(A) fluxes of 
the emission lines, and the extinction coefficients C(H/3) are 
shown in TableQ] The extinction coefficients C(H/3) = 0.74, 0.27, 
and 0.31 for Haro 11B, C, and ESO 338-IG 004, respectively, 
are derived from the hydrogen Balmer lines only in the UVB 
and VIS arm spectra. They are adopted for the correction of 
emission-line fluxes not only in the UVB and VIS spectra, but 
also in the NIR range. We also show in Table[T]the theoretical re- 
combination fluxes for hydrogen Balmer, Pa schen, and Brackett 
lines calculated by Hummer & Storey (1987) (last column of the 
Table), adopting a case B model with an electron temperature T e 
= 12500 K, and an electron number density N e = 100 cm~ 3 . 

The deviations of the observed fluxes of bright Balmer, 
Paschen and Brackett hydrogen lines after correction for extinc- 
tion with a single value of C(H/3) from the theoretical recom- 
bination values (Hummer & Storev, 1987) do not exceed a few 
percent (cf. Table [TJ. The exceptions are extinction-corrected 
fluxes of several Paschen and Brackett hydrogen lines marked in 
Table [TJ which deviate significantly from the theoretical values 
(e.g., PayS line) because of strong telluric absorption. Excluding 
these lines, we may conclude that a single C(H/3) value derived 
from the hydrogen Balmer lines can be used over the entire 
wavelength range of -AA3000 - 24000. 

Thus, we confirm our previous findings obtained for sev- 
eral low-metallicity emission-line galaxies (Mrk 59, II Zw 40, 
Mrk 71, Mrk 996, SBS 0335-052E, PHL 293B, and GRB 
HG 031203) that the extinction coefficient C(H/3) is not sys- 
tematically higher across the NIR wavelength range than the 
optical range and that no additional regions contribute to 
the NIR l ine emission that are, however, hidden in the opti- 
cal range dVanzi et all l200d l2002t lizotov et all [20091 1201 lbt 
llzotov & ThuanLl201 UlGuseva et alll201 lbl) . 



of log N/O, log Ne/O, and log S/O for ESO 338-IG 04 (Tabled 
are consistent with the data of Masegosa et al. ( 1994) (log N/O = 
-1.40±0.01 and log Ng/D = -O75+0. 10). lRaimann et all J2000) 
(log N/O = -1.57), Bergvall (19 81 (lo g N/O = -1.18, and log 
S/O = -1.63) lBergvall& Ostlin (2002) (log N/O = -1.5), and 
iGuseva et all (1201 lal) (log N/O = -1.34, log Ne/O = -0.78, and 
log S/O = -1.60). 

3.2. Hidden star formation 

The spectra of the studied galaxies were obtained simultaneously 
over the entire optical and near-infrared wavelength ranges. This 



3.3. H 2 ro-vibrational emission 

Molecular hydrogen lines do not originate in H II regions, but 
in neutral molecular clouds. In the NIR, they are excited by 
two mechanisms. The first is a thermal mechanism consisting of 
collisions between neutral species (e.g., H, H2), resulting from 
large-scale shocks driven by powerful stellar winds, supernova 
(SN) remnants, or molecular cloud-cloud collisions. The second 
is fluorescence due to the absorption of ultraviolet photons pro- 
duced by hot stars. By comparing the obser ved line ratios with 
those predicted by for instance the models of Black & Dishoeck 
d!987l) . it is possible to discriminate between the two processes. 
In particular, line emission from the vibrational level v = 2 and 
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Fig. 6. Best-fit model SED fits to the redshift- and extinction- 
corrected observed spectra of Haro 11B, Haro 11C, and ESO 
338-IG 004. The contributions from the stellar and ionised gas 
components are shown by the green and blue lines, respectively. 
The sum of stellar and ionised gas emission is shown by the red 
line. 



higher vibrational levels are virtually absent in collisionally ex- 
cited spectra, while they are relatively strong in fluorescent spec- 
tra. 

Thirteen ro-vibrational H2 emission lines are detected in the 
NIR spectrum of Haro 11B (see Fig. [TJ. Only two H2 emis- 
sion lines, 1.957/zm 1-0 5(3) and 2.122/mi 1-0 5(1), are de- 
tected in the spectrum of ESO 338-IG 004 and no H2 lines 
are detected in the spectrum of Haro 11C (Table [T). The fluxes 
of the H2 lines relative to that of the strongest 2.122 /im 1- 
5(1) line are shown in Table |4] Two H2 emission lines, the 



1.188 /mi 2-0 5(0) and the 1.238 /mi 2-0 g(l) transitions are 
blended with the Fe II lines. The very large observed flux ra- 
tio of the 1.188 fj.m line (Table |4| can be due to contamina- 
tion by a singly ionised iron line. T he last two column s show 
the theoretical ratios calculated by Black & Dishoeck (1987) 
in the cases of fluorescent and collisional excitation. It can be 
seen that the observed line flux ratios for most of the H2 lines 
agree with those predicted for fluorescent excitation, support- 
ing the previous findings that fluorescence is the main excitation 
mechanism of ro-vibrati onal H2 fine s in most BCDs with high- 
excitat i on H 11 regions dVanzi et all l2000l l2008t llzotovetaH 
2009a; Izotov & Thuan, 1201 lb. However, t here are still some 
BCDs (e. g., PHL 293B llzotovet alll2011bh for which the col- 
lisional excitation of ro-vibrational H2 lines is the dominat exci- 
tation. 

3.4. CLOUDY models 



Using the CLOUDY code (version c 10.00) of iFerland et all 
( 1998), we examine the excitation mechanisms not only of the 
hydrogen lines but all emission lines by producing stellar pho- 
toionisation models for our objects. The CLOUDY input param- 
eters for the best-fit models with a single stellar population are 
shown in Table [3] The logarithm of the number of ionising pho- 
tons per second Q(H) is calculated from the extinction-corrected 
Hy6 luminosity. We adopt the shape of the ionising stellar radi- 
ation of the Starburst-99 models (Leithereret al., 1999) with a 
starburst age shown in Table [5] The electron number density N e 
is assumed to be constant with radius. The ratios of abundances 
of different species to hydrogen are taken from this paper, except 
for carbon which has no bright emission lines in the optical and 
NIR ranges. To obtain the abundance of this species, we use a 
relation C/O v s. oxygen abundanc e for the emission-line galax- 
ies studied bv lGarnett et all dl995l) . The emission-line fluxes for 
the best-fit models are shown in Table [6] 

The CLOUDY-predicted and extinction-corrected fluxes of 
bright emission lines in Haro 11B and Haro 11C are shown 
in Table [6] The overall agreement achieved implies that a H 
II region model including only the stellar ionising radiation of 
a single stellar population is able to account for the observed 
fluxes, in both the optical and NIR ranges. No additional ex- 
citation mechanism such as shocks from stellar winds and su- 
pernova remnants is needed. This conclusion also holds for the 
[Fe II] A 12570, 16436 emission lines ( Table |6]), which are often 
considered as shock ind icators (e.g. iMoorwood & Oli va. 1988; 
IOHvaetalllT990Ll200lh . 

We note that the photoionisation models for a single stellar 
population cannot reproduce the extinction-corrected fluxes in 
ESO 338-IG 004, because they are measured in the spectrum 
that includes two star-forming regions with different excitation 
properties. The superposition of two CLOUDY models should 
be used instead to fit the observations of ESO 338-IG 004. 

To extend our study of Ha ro 11 into the M IR range, we 
used the Spitzer observations of IWu et all ([2008) obtained with 
an aperture, which is much larger than the one used in the X- 
shooter observations, and includes all of the bright H II regions 
A, B, and C. Therefore, to compare the observed and CLOUDY- 
predicted emission-line fluxes in the MIR range (relative to H/3) 
we should take into account aperture corrections. In particular, 
the extinction-corrected H/3 flux from the X-shooter observations 
cannot be used because it was obtained with a smaller aper- 
ture. We instead derive the H/3 extinction-corrected flux from 
th e total Ha flux F (Hg)=3.326xl0~ 12 erg s -1 crrT 2 measured 
bv lWu et all ([2008) after its correction for extinction. We do not 
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Table 4. H2 emission line flux ratios for Haro 1 IB 



,1 
yum 



Line 



Observations" 



Model" 



Fluor. Collis. 



1.162 
1.185 
1.188 
1.207 
1.233 
1.238 
1.601 
1.957 
2.034 
2.073 
2.122 
2.223 
2.248 



H 2 2-l 
H 2 3-l 
H 2 2-0 
H 2 3-l 
H 2 3-l 
H 2 2-0 
H 2 6-4 
H 2 1-0 
H 2 1-0 
H 2 2-l 
H 2 1-0 
H 2 1-0 
H 2 2-l 



5(1) 

5(3) 

5(0)+[FeIl] 

5(2) 

5(1) 

fi(l)+Fe II? 

2(1) 

5(3) 
5(2) 
5(3) 
5(1) 
5(0) 
5(1) 



0.2 
0.2 
0.9 
0.1 
0.6 
0.2 
0.1 
0.5 
0.3 
0.3 
1.0 
0.2 
0.2 



0.4 
0.4 C 
0.2 C 

0.2 
0.5 
0.4 
0.4 
0.7 C 
0.5 
0.2 
1.0 
0.6 
0.5 



0.0 



0.0 
0.0 
0.0 
0.0 

0.3 
0.0 
1.0 
0.3 
0.0 



"Flux ratios relative to the H 2 2. 122 jum flux. 

b Model values by Black & Dishoeck 1 1987). We adopt their model 1 
for fluorescent lines and model SI for collisionally exc i ted lin es. 
c Flux ratios are from model 14 by Black & Dishoeck ( 119871) . because 
these lines are not present in model 1 . 



Table 5. Input parameters for the CLOUDY model 



Parameter 


Haro 1 IB 


Haro 11C 


log Q(Hf 


53.51 


52.30 


Starburst age, log Myr 


6.35 


6.50 


Y e , cirT 3 


10 


10 


f 


0.02 


0.06 


log He/H 


-0.93 


-0.93 


log C/H 


-3.97 


-4.20 


log N/H 


^1.58 


^1.69 


logO/H 


-3.67 


-3.90 


log Ne/H 


^.33 


^.53 


log S/H 


-5.50 


-5.64 


log Ar/H 


-6.17 


-6.22 



"2(H) is the number of ionising photons per second. 
b Volume filling factor. 



know the average extinction coefficient C(H/3) within an aper- 
ture used for the MIR observations by IWu et al.l ([2008). It is 
known from HST data (Ad amoet all 1201 Ol) that the redden- 
ing in Haro 11 varies in the range E(B - V)=0.0 to 0.6 (cor- 
responding to C(HjS)=0.0 to 0.88), and to even higher values, 
in agreement with our extinction determinations for Haro 1 1 B 
and Haro 11C. Therefore, to correct the H/3 flux for extinction 
within a large aperture, we consider two limiting cases by adopt- 
ing the C(Hj6) values derived from both the Haro 1 IB spectrum 
(high extinction) and the Haro 11C spectrum (low extinction). 
It is also reasonable to assume that the averaged extinction in 
the area covered by the MIR observations is between these two 
limiting cases. In the third and fourth columns of Table [6] we 
show the MIR emission-line fluxes normalised to the extinction- 
corrected HjS flux with C(HjS)=0.745 (the case of Haro 1 IB) 
and C(H/3)=0.27 (the case of Haro 11C), respectively. We also 
note that a correction of MIR emission lines for extinction is not 
needed. In the last column of Table [6] for MIR data, we show 
CLOUDY predictions with input parameters from Table |5] It is 
seen from Table [6] that the predicted relative MIR emission-line 
fluxes are between the two limiting cases of extinction-corrected 
fluxes with high and low extinctions, suggesting that with a more 
realistic intermediate extinction the agreement between models 
and observations would be much better. This implies that the 



Table 6. Comparison of the observed and predicted intensities 
of some strong emission lines in Haro 1 IB and C 



Ion 



CLOUDY 



Observations" 



B* 



£5 



B* 



C^ 



///(H/3) ///(H/3) ///(H/3) ///(Hfl) 



3727 [O II] C 


2.058 


1.927 


2.186 


2.336 


3869 [Ne ill] 


0.283 


0.321 


0.299 


0.314 


4101 H5 


0.264 


0.267 


0.277 


0.269 


4340 Hy 


0.472 


0.475 


0.499 


0.461 


4861 H0 


1.000 


1.000 


1.000 


1.000 


4959 [O III] 


1.208 


1.070 


1.194 


0.990 


5007 [O in] 


3.635 


3.220 


3.716 


3.172 


6563 Ho- 


2.917 


2.891 


2.907 


2.864 


6583 [N 11] 


0.342 


0.381 


0.488 


0.576 


6716 [S 11] 


0.072 


0.078 


0.156 


0.291 


6731 [S 11] 


0.051 


0.055 


0.102 


0.134 


9069 [S in] 


0.172 


0.154 


0.156 


0.144 


9532 [S in] 


0.426 


0.382 


0.432 


0.385 


10829 He 1 


0.286 


0.254 


0.326 


0.283 


12570 [Fe II] 


0.021 


0.030 


0.014 


0.028 


12821 Pa£ 


0.164 


0.161 


0.144 


0.176 


16436 [Fe II] 


0.018 


0.025 


0.012 


0.017 


19450 Bx6 


0.018 


0.018 


0.019 


0.014 


21661 Bry 


0.028 


0.027 


0.031 


0.028 



MIR data for Haro 1 1 


Ion 


F A 


///(H/3) c 


///(H/3)' 


CLOUDY 


10.51;um[SlV] 
12.81;um[NeIl] 
15.55^m [Ne III] 
18.71yum [S III] 


4.395 
3.123 
9.814 
4.524 


0.122 
0.087 
0.273 
0.126 


0.365 
0.259 
0.814 
0.376 


0.224 
0.097 
0.466 
0.206 



"extinction-corrected fluxes. 

b B and C mean Haro 1 IB and Haro 1 1C, respectively. 

c sum of the 3726 and 3728 [O II] lines. 

d observed flux in units of 10" 13 erg s -1 

c extinction-corrected flux ratios with C(Hy6)=0.745. 

'extinction-corrected flux ratios with C(HjS)=0.27. 



cm" 2 JWuetall l2008). 



MIR range, similar to the NIR range, does not resemble highly 
extincted emission-line regions that significantly contribute to 
the MIR line emission. 



3.5. Evidence for luminous blue variable (LBV) stars in 
Haro 1 1 

The Ha emission line in the spectrum of Haro 11C is asym- 

1 ■■ 1 I I 

metric and has also a broad component. Ostlin et al. (2001) also 
reported that the Ha line profiles of Haro 1 1 are broad, up to full 
widths at half maximum (FWHM) of 270 km s ' , and have a non- 
Gaussian shape. They suggested that two or more non-virialised 
ionised gas components may be present in the galaxy. The 
extinction-corrected relative Ha flux /(Ha)//(Hy6) = 3.82 in Haro 
11C, obtained with C(H/?)=0.27, deviates significantly from 
the th eoretical recombination value 2.86 dHummer & Storevl 
dl987D for the electron temperature 10000 K and the electron 
number density N e = 100 cirT 3 when the observed Ha flux 
F(Ha)=3.350xlO~ 14 erg s~' cirT 2 is adopted. This high Ha flux 
may be due to substantial emission from either supernovae of 
type II or circumstellar envelopes of stars. Because of its anoma- 
lous behavior this line was excluded in sect. 13. II from the cal- 
culation of the extinction coefficient. The decomposition of the 
Ha emission line of Haro 11C into three Gaussian emission 
components and one Gaussian absorption component is shown 
in Fig. [5] Two of the Gaussian emission components (the nar- 
rowest and the broadest) marked by blue dotted lines are likely 
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Table 7. Eqiuvalent widths of the stellar absorption CO band 
(A=23fj.m) in Haro 1 1 



sent in the s pectra of the Orion neb- 
etal.1 (l2000h . lEsteban et all (|2004), and 



Starburst-99 a 



inst. c 



cont. 



Observations 
Haro 1 IB 8 
Haro 11C 9 



Z=0.004 Z=0.008 Z=0.004 Z=0.008 



8-11 



13-15 



8-10 



"Values for the time of 8-10 Myr after the onset of star formation 

dLeithereretal.1119991) . 

b Z is heavy element mass fraction. 

instantaneous burst model. 

d Continuous star formation. 



linked to the nebular emission, implying that there is low- and 
high-velocity ionised gas in the H II region. These two com- 
ponents are sufficient to reproduce the shape of the strongest 
forbidden line [O III] /15007, strengthening the evidence for a 
nebular origin of these components. However, the Ha profile is 
more complex, and two additional components (emission and 
absorption) are needed. The green dotted line is the composi- 
tion of the third emission and absorption components and rep- 
resents a P Cyg profile. The observed Ha profile (black solid 
line) is nicely fit by the summed profile, consisting of two neb- 
ular components and the P Cyg profile (red dashed line). The 
flux F(Ha)=2.281xl(T 14 erg s" 1 cirT 2 of the Ha line after sub- 
traction of the P Cyg profile corresponds to the recombination 
ratio /(Ha)//(Hj6) = 2.86 (Table [B after the correction for ex- 
tinction with C(H/3) = 0.27 is applied. Thus, the inclusion of the 
stellar component of Ha emission is required to fit the observa- 
tions, implying that either LBV star(s) or supernovae are present. 
The Ha emission line is therefore the most affected by the pres- 
ence of a stellar wind. The observed Balmer decrement of the 
broad emission lines in LB Vs and SNe is often much larger than 
that for the narrow lines, reflecting the high density and possibly 
high extinction in the circumstellar envelopes where the broad 
lines originate. In particular, this is true for the LBV in PHL 
293B and resu lts in a very strong Ha emission line comp ared 
to the H/3 line dlzotov & Thuanl 120091: Hzotov et all 1201 lbl) . On 
the other hand, the small FWHM of ~ 180 km s" 1 of the P Cyg 
component precludes the supernovae scenario. We note that, in 
principle, in ESO 338-IG 004 some additional contribution from 
LBV stars or shocks to the Ha emission can also be present. 
Therefore, C(H/3) was derived excluding Ha, hence /(Ha)//(H/3) 
= 2.94 (Table[TJi is slightly higher than the recombination ratio. 

The extinction-corrected luminosity of the Ha P Cyg com- 
ponent in Haro 11C is very high, L(P Cyg) = 1.32xl0 40 erg 
s _1 , exceeding by a factor of ~ 10 - 100 the Ha luminosi- 
ties of LBV stars in the low-metallicity emission-line galaxies 
NGC 2366 (Drisse n et al.lll997ll200lh, PHL 293B, and DDO 6 8 
dlzotov et aUl2007l:IIzotov & Thuanl 120091: llzotov et alll201 lbl) . 

Moreover, sixteen permitted Fe II emission lines are ob- 
served in the spectra of Haro 11B and Haro 11C (see Figs. |4] 
[U |2] [3] and Table [TJ>, which reflects their circumstellar origin 
in the dense stellar wind. No such lines are present in low- 
density nebular gas. The Orion Nebula serves as a template for 
the investigation of different excitation mechanisms of emis- 
sion lines observed in a gaseous nebulae. It is one of the most 
well-studied of nearby Galactic H II regions with an extensive 
list of ide ntifications of various lines, i nclud i ng the very weak 
ones (e.g. Peimbert & Torres -Peimbertl Il977t lOsterbrock et all 
1992). There are high-resolution spectroscopic observations of 
the Orion nebula with high signal-to-noise ratios. Nevertheless, 
the Fe II lines observed in the VLT/X-shooter spectra of Haro 



11B and C ar e not present in the s pectra of the Orion neb- 
ula obtained b y Baldwin e 
iMesa-Delgado et alj d20~ 

On the other hand, the permitted Fe II emission lines are 
common in spectra of LBV stars undergoing giant eruption 
events during the late evolutionary stages of the most massive 
stars. One of the most well-known LBVs is the binary sys- 
tem 77 Carina (77 Car), which undergoes periodical outbursts 
(e.g. [M assey, 2003; Smith et al., 2004; Damineli et al., 2008 
IClark et aUl2005H Damineri et al., 2008at lSoker & Kashill201 
The bright emission lines of H, He, and the permitted Fe II and 
forbidden [Fe II] lines are observed in LBVs that also display 
P Cyg profiles during eruption events. Most of the emission 
lines of different species (Mg I, [Fe II], Si II, [Ni II], [Cr II], 
N I, and [C I]) (Table [TJ are observed in nebulae as well as in 
emission-line stars like LBVs, supergiants with H and He emis- 
sion, and Of/WN9 stars. We also note that all Fe II emission 
lines seen in the spectra of Haro 1 IB and Haro 1 1 C are detected 
only in the hig h-resolution sp ectra of 77 Car by e.g. lHamannet alj 
d!994h and lWallerstein et al.1 (120011). For the line ide ntification, 
we use the data for 77 Car of Wallerstein et al. (2001) in a wide 
wavelength range ~/13100 - 9000 (with several gaps due to ob- 
servational program reasons or telluric absorption) and data of 
Ham ann etall (11994 in the wavelength range -/16450 - 24500. 
These two sets of data cover the whole wavelength range of the 
VLT/X-shooter observations excluding gaps. Additiona lly, fo r 
some line identifi cati ons we used the d ata of iGaviolal d!953l) . 
Uohanssonl d!977l) . and lAllen et ail d!985l) . 

We find that LBV is probably also present in Haro 1 IB. This 
is indicated by the Fe II and [Fe II] emission lines in its spectrum 
(Fig. |U Table [TJ). It is also possible that broad Ha emission is 
also present, but masked by very strong nebular Ha and [N II] 
/16548, 6583 emission lines. On the other hand, no permitted Fe 
II emission lines were detected in the spectrum of ESO 338-IG 
004 (see TableQ]). 

Thus, we have three pieces of evidence supporting the idea 
that a LBV star is present in Haro 11C. These are that a) a P 
Cyg profile is needed to reproduce the observed Ha profile, b) 
the only possible way to adjust the observed Balmer decrement 
to the theoretical one is the subtraction of the P Cyg component 
from the Ha emission line, and c) the presence of numerous Fe 
II emission lines, which is a common feature of the LBV star 
spect rum. 

I Smith et aD d201lh collected data for a number of extragalac- 
tic optical transients or "supernova impostors" related to giant 
eruptions of LBVs and illustrated their remarkably wide diver- 
sity in terms of their peak absolute magnitude, duration, pro- 
genitor stars, outburst spectra, and other observable properties. 
In particular, the peak absolute magnitude varies from -10 mag 
to -16 mag and increases during eruption by more than 2-4 
mag for LBV stars of the Milky Way (77 Car, P Cyg), VI in NGC 
2366, HD 5980 in the SMC, and others. There is not much infor- 
mation in the literature about the range of LBV Ha luminosities, 
but it is likely that they have a large range, which is similar to the 
absolute magnitudes. Therefore, it is difficult to quantitatively 
estimate the number of LBV stars in Haro 1 1C. 



3.6. Spectral energy distributions, stellar masses and a 
red continuum excess 

We derive the stellar masses of three H II regions (see Table [3) 
by fitting their spectra to a superposition of SEDs of single stel- 
lar populations with differe nt a ges following the pres criptions of 
iGuseva et alJ (2006, 20071) and llzotov etall d201 lal) . This tech- 
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nique also allows us to check whether a red continuum excess 
is present in the spectra of Haro 11B and Haro 11C. We note, 
that for star- forming galaxies with EW(H/3) > 100 the ionised 
gas continuum emission contributes appreciably to the total con- 
tinuum emission. Neglecting it can lead to an overestimation of 
galactic stellar masses by a factor of a few ( Izotov et al ., 201 la). 
We include the ionised gas emission in the fit of the spectra of 
all three objects, despite the equivalent width of the H/3 emission 
line being large (EW(H/3)=155) only in the Haro 1 IB spectrum, 
while EW(Hy6) < 100 in the spectra of Haro 1 1C and ESO 338- 
IG 004 (TableQ). 

As each SED is the sum of both stellar and ionised gas emis- 
sion, its shape depends on the relative contributions of these two 
components. In active star-forming regions, the contribution of 
the ionised gas emission can be very large. However, the equiv- 
alent widths of the H/3 emission line never reach the theoreti- 
cal values for pure gaseous emission of ~ 900 - 1100, imply- 
ing that there is a non-zero contribution of stellar emission in 
all objects. To take into account gaseous emission, we use the 
observed emission-line fluxes and equivalent widths. The con- 
tribution of the gaseous emission is scaled to the stellar emission 
using the ratio of the observed equivalent width of the H/3 emis- 
sion line to the equivalent width of H/3 expected for pure gaseous 
emissi on. The continuous gaseous emission is taken from Aller 
(1984) and includes hydrogen and helium free-bound, free-free 
and two-photon emission. In our models, this emission is always 
calculated using the electron temperature T e (H + ) of the H + zone 
and with the He/H abundance ratio derived from the H II region 
spectrum. 

To calculate the contribution of stellar emission to the 
S EDs, we adop t a grid of the Padua stellar evolution models 
of lGirardiet al. (2000|3 with a heavy element mass fraction Z 
= 0.008 for Haro 11B and Z = 0.004 for Haro 11C and ESO 
338-IG 004. Using these data, we calculate w ith the package 
PEGASE.2 dFioc & Rocca-Volmerangel 1 19971) a grid of stellar 
instantaneous burst SEDs for a range of ages from 0.5 Myr to 
15 Gyr. We adopt a stellar initial mass function with a Salpeter 
slope, an upper mass limit of 100 M , and a lower mass limit of 
0.1 M . The SED for any star formation history (SFH) can then 
be obtained by integrating the instantaneous burst SEDs over 
time with a specified time-varying star formation rate. We ap- 
proximate the SFH by a combination of a recent short burst to 
account for the young stellar population, and continuous star for- 
mation to describe the properties of the older stars. 

Finally, the contribution of the nebular continuum and emis- 
sion lines were taken into account. Each fit was performed over 
the whole observed spectral range ~AA3QQQ - 24000. We per- 
formed 5xl0 5 Monte Carlo models by varying simultaneously 
the age of young stellar population f(young), the age of the old- 
est stars f(old), the mass ratio of old-to-young stellar populations 
b = M(old)/M(young), the electron temperature r e (H + ) in the 
H + zone, and the extinction coefficients for gas and for stars, 
C(HyS) gas and C(HjS) stars , respectively. 

The best-fit model SED is found from^ 2 minimisation of the 
deviation between the modelled and the observed continuum in 
ten ranges of the whole spectrum free of emission lines. 

In Fig.|6l we show the best-fit model SED fits to the redshift- 
and extinction-corrected spectra of Haro 11B, C, and ESO 338- 
IG 004. The contributions from the stellar and ionised gas com- 
ponents are shown by green and blue lines, respectively. The sum 
of both stellar and ionised gas emission is shown by a red line. 
It is seen that the contribution of gaseous emission is essential 
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only for Haro 11B because of the high EW(HyS). For an other 
two objects, this contribution is small, especially for Haro 11C 
wit h EW(Hff) = 1 6. 

lAdamo et al.l (l2010h found using HST photometry data for 
Haro 11, that a strong red excess is present in spectra of Haro 
11B and Haro 11C at wavelengths > 8000 with respect to the 
synthetic evolutionary models. They show that despite the tighter 
model fit to the Haro 1 IB spectrum achieved (y 2 for Haro 1 IB 
is smaller than that for Haro 11C), the red excess in Haro 11B 
is larger. Fig. [6] shows that all observed spectra in the whole 
range of wavelengths ~/M3000-24000 are fitted quite well by 
the model SEDs, except for the spectrum of Haro 11B where 
the flu x excess at wavelengths >1.6/im is clearly present. Thus, 
Adamo et al. (2010) found red excesses in the B and C knots but 
in the present paper we find one only for the B knot. This dis- 
agreement can be partially associated with the different sizes of 
the extracted regions (we ad opted extraction aper tures of ~l"x 
4" for Haro 1 IB and C, while Adamo et al. (2010) adopted much 
smaller apertures). However, we note that the majority of the 
light in the B and C knots within the large adopted spectroscopic 
apertures comes from the brightest compact clusters within the 
slit, implying that we measure simi lar continuum flux d ensities 
to those in the photometric study of Ada mo et al.l (2010). 

What is the nature of this excess? As mentioned above the 
observed fluxes of the Balmer, Paschen, and Brackett hydrogen 
lines after correction for extinction with a single value C(HyS) 
(see Table [TJ agree well with the theoretical recombination val- 
ues of Hummer & Storev (1987) for all studied star-forming re- 
gions, i.e. Haro 11B, Haro 11C, and ESO 338-IG 004. The 
same result was obtained from CLOUDY modelling (Table [6}. 
This implies that the NIR excess is not due to highly obscured 
hot stars ionising the interstellar medium, which would require 
that the continuum excess be followed by an excess of the NIR 
emission lines, a result we do not find. Evidently, only the non- 
ionising cool stars, including a large population of luminous red 
supergiant stars (RSGs), can contribute to the NIR continuum 
excess. 

This conclusion is supported by the detection of the strongest 
stellar absorption CO band at ~2.3 fim in the spectra of Haro 1 IB 
and Haro 11C. This spectral feature is a signature of red giant 
and red supergiant stars. The equivalent width of the CO band 
depends on stellar temperature and luminosity and decreases 
with both increasing stellar temperature and decreasing stellar 
luminosity. We measure equivalent widths EW(CO 2.3 fim) of 8 
and 9 in Haro 11B and Haro 11C, respectively. In Table [7] we 
compare the observed value s with the Starburst-99 stellar pop- 
ulation synthesis models dLeitherer et all Q999) and conclude 
that the stellar population of an older burst with an age of ~ 10 
Myr including red supergiant stars can be responsible for the ob- 
served red excess and CO absorption at ~2.3 fim. This older non- 
ionising stellar population is probably more obscured by the dust 
and therefore does not contribute to the optical light. However, 
we are unable to infer whether this is true from the hydrogen- 
line decrement because the older population does not produce 
emission in hydrogen lines. 

We derive the instantaneous S FR(Ha), using the extinction- 
corrected luminosity L(H a) of the Hg emis sion line and the fre- 
quently used relation by iKennicutn (JT998), which, however, is 
obtained for longer-duration star formation in disc galaxies and 
may not be valid for knots B and C where stars were formed 
almost instantaneously. We obtain S FR = 3.5, 0.22, and 0.11 
M Q yr 1 in Haro 1 IB, Haro 1 1C, and ESO 338-IG 004, respec- 
tively, for the regions within the apertures used for the spectro- 
scopic observations (see sect.|2]i. To achieve this result, we adopt 
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the distances of 82.4 Mpc and 38.4 Mpc for Haro 1 1 and ESO 
338-IG 004, respectively, which are taken from the NED. The 
star formation rates o f luminous compact emission-line galax- 
ies (LCGs) derived bv llzotovet al.l (1201 lal) are in the range 0.7 
- 60 M Q yr~' with an average value of ~4 M© yr , Thus, Haro 
1 1 is among the galaxies with significant star-formation activity 
in its H II region B. The specific star formation rates {SSFR) 
are then defined as S S FR(Ha) = S FR(HayM„ and are equal to 
2.8xl0~ 8 yr" 1 , l.lxlO -9 yr 1 , and 2.0xl0~ 9 yr 1 for Haro 11B, 
Haro 11C, and ESO 338-IG 004, respectively. In addition, the 
S SFR for Haro 11B is high and similar to those found in LCGs 
by llzotov et al.l d2011al) and in high-redshift galaxies, while the 
S SFRs for Haro 1 1C and ESO 338-IG 004 are lower by a factor 
of- 10. 

4. Conclusions 

We have studied the spectra of the blue compact galaxies Haro 
1 1 and ESO 338-IG 004 obtained by performing VLT/X-shooter 
spectroscopic observations in the wavelength range ~AA3000 - 
24000. We have arrived at the following conclusions: 

1. We derived oxygen abundances of 12 + log O/H = 8.33 + 
0.01 and 8.10 + 0.04 in the two H II regions Haro 1 IB and Haro 
11C, respectively, w hich are appreciably high er than the previ- 
ous determination of iBergvall & Ostlinl d2002l) (12 + log O/H = 
7.9). Our oxygen abundance 12 + log O/H = 7.89 ± 0.01 for 
ESO 338-IG 004 is slightly lower than t he previous determina- 
tions, which are in the range 7.9 2 - 8.08 (Masegosa et al., 1994; 
BerevallL 119851: iRaimann et all l2000t ^Bergvall & OsflinL 12002: 
Gusevaet all 1201 lal) . 

2. We used X-shooter data together with Spitzer observations 
in the mid-infrared range to test for hidden star formation. The 
observed fluxes of hydrogen lines correspond to the theoretical 
recombination values after correction for extinction with a single 
value of C(HjS) applied to a wide wavelength range from near- 
UV to NIR and MIR. Thus, we confirm our previous findings 
obtained for several low-metallicity emission-line galaxies (Mrk 
59, II Zw 40, Mrk 71, Mrk 996, SBS 0335-052E, PHL 293B, 
and GRB HG 031203) that the extinction coefficient C(B/3) is 
not systematically higher in the NIR wavelength range than the 
optical one and that no hidden star formation contributes appre- 
ciably to the NIR line emission but not the optical line emission. 

3. All observed spectra are fitted quite well across the entire 
wavelength range covered by model SEDs, except for the con- 
tinuum flux excess at wavelengths >1.6/im for Haro 11B. An 
anomalously large population of luminous red supergiant stars 
(RSGs) in Haro 11B is probably responsible for the observed 
flux excess relative to the model SED in the NIR range. 

4. Thirteen ro-vibrational H2 emission lines are detected in 
the NIR spectrum of Haro 11B. From the observed flux ratios 
of the H2 lines, we conclude that fluorescence is the main exci- 
tation mechanism of ro-vibrational H2 lines. This confirms the 
previous finding for BCGs with high-excitation H II regions. 

5. The agreement between the CLOUDY-predicted and 
extinction-corrected emission-line fluxes implies that a H II re- 
gion model including only stellar photoionisation as an ionising 
source is able to account for the observed fluxes in the optical, 
NIR, and MIR ranges. No additional excitation mechanism such 
as shocks from stellar winds and supernova remnants is needed. 

6. We find evidence of an LBV star in Haro 11C: a) a P 
Cyg profile is needed to reproduce the asymmetric Ho- emission 
line; b) the same P Cyg emission is needed to explain a strong 
Ha emission excess above the recombination value; c) numer- 
ous permitted singly ionised iron emission lines are detected in 



the spectrum of Haro 1 1C. These lines are common in the spec- 
tra of LBV stars. The high luminosity of the P Cyg He line L(P 
Cyg) = 1.3xl0 40 erg s _1 exceeds by a factor of 10 - 100 the 
LBV luminosities in some low-metallicity emission-line galax- 
ies. Similarly, the numerous Fe II and [Fe II] emission lines in 
the spectrum of Haro 1 IB indicates that LBV is also present in 
this knot. We find no evidence of LBV in ESO 338-IG 004. 
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Fig. 3. Same as Fig.Q]but for ESO 338-IG 004. 
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Table 1. Emission Line Intensities 









Galaxy 










Haro 1 1 B 


Haro 11C 


ESO 338-IG 004 
F(i)/F(H^) a 7(i)//(H/3) D 


theory 


Ion 


F(A)/F(Hfif 


KWiW? 


F(i)/F(H^) a 


I{X)II(Wf 


Case B c 








a) Near-UV and optical range 








3188 He I 


2.14 ±0.07 


4.92 ±0.16 


3.01 ± 0.32 


3.99 ± 0.43 








3322 [Fe III]? 


0.21 ± 0.04 


0.45 ±0.10 












3614 He I 


0.19 ±0.03 


0.34 ± 0.05 












3634 He I 


0.15 ±0.03 


0.26 ± 0.05 












3671 H24 


0.16 ±0.04 


0.31 ±0.08 






0.19 ±0.04 


0.24 ± 0.05 


0.39 


3676 H22 


0.24 ± 0.04 


0.40 ± 0.07 






0.42 ± 0.05 


0.53 ± 0.06 


0.50 


3679 H21 


0.25 ± 0.03 


0.43 ± 0.05 






0.34 ± 0.04 


0.43 ± 0.06 


0.57 


3683 H20 


0.30 ± 0.02 


0.51 ±0.03 










0.66 


3687 H19 


0.34 ± 0.02 


0.59 ± 0.03 










0.77 


3692 H18 


0.49 ± 0.02 


0.85 ± 0.03 










0.91 


3697 H17 


0.61 ± 0.02 


1.05 ±0.04 










1.08 


3704 H16 


1.06 ±0.03 


1.82 ±0.05 










1.29 


3712 H15 


0.97 ± 0.03 


1.66 ±0.05 










1.57 


3722 H14 


1.82 ±0.04 


3.09 ± 0.07 


1.41 ±0.14 


1.67 ±0.17 






1.93 


3726 [0 II] 


63.69 ± 0.91 


107.80 ± 1.68 


85.01 ± 1.25 


100.74 ± 1.64 


46.53 ± 0.70 


57.82 ± 0.94 




3728 [O II] 


65.57 ± 0.93 


110.83 ± 1.73 


112.19 ± 1.65 


132.86 ±2.16 


64.12 ±0.95 


79.64 ± 1.28 




3733 He I 


2.21 ±0.13 


3.72 ± 0.23 












3734 H13 


1.62 ±0.10 


2.73 ±0.17 


1.41 ±0.27 


1.66 ±0.33 


1.39 ±0.07 


1.72 ±0.09 


2.41 


3750 H12 


1.72 ±0.05 


3.35 ±0.13 


1.56 ±0.24 


4.95 ± 0.99 


2.28 ±0.12 


3.35 ± 0.24 


3.07 


3771 HI 1 


2.20 ± 0.06 


4.12 ±0.14 


2.29 ± 0.34 


5.87 ± 1.06 


2.66 ±0.12 


3.80 ± 0.24 


4.00 


3798 H10 


3.11 ±0.05 


5.56 ±0.13 


3.18 ±0.17 


6.87 ± 0.52 


3.92 ±0.18 


5.29 ± 0.29 


5.34 


3820 He I 


0.56 ± 0.02 


0.90 ± 0.03 


0.11 ±0.08 


0.12 ±0.10 








3835 H9 


4.66 ± 0.07 


7.91 ±0.15 


4.31 ±0.15 


8.03 ± 0.44 


5.84 ±0.13 


7.57 ± 0.23 


7.37 


3869 [Ne III] 


19.07 ± 0.27 


29.96 ± 0.46 


27.21 ± 0.45 


31.38 ±0.56 


39.79 ± 0.59 


47.93 ± 0.76 




3889 He I + H8 


12.07 ±0.17 


19.22 ± 0.31 


11.65 ±0.24 


16.54 ± 0.48 


16.52 ± 0.26 


20.29 ± 0.37 


10.60 


3967 [Ne III] 


6.08 ±0.17 


9.08 ± 0.26 


6.25 ± 0.29 


7.08 ± 0.34 


11.05 ±0.17 


13.04 ± 0.21 




3970 H7 


10.11 ±0.16 


15.10 ±0.25 


11.20 ±0.34 


12.69 ± 0.40 


12.50 ± 0.22 


14.74 ± 0.27 


16.00 


4026 He I 


1.10 ±0.02 


1.60 ±0.04 


0.74 ±0.17 


0.82 ± 0.20 


1.20 ±0.10 


1.40 ±0.11 




4068 [S II] 


1.52 ±0.03 


2.15 ±0.05 


5.15 ±0.26 


5.73 ± 0.30 


0.84 ± 0.06 


0.98 ± 0.07 




4076 [S II] 


0.45 ± 0.02 


0.64 ± 0.03 












4101 H5 


19.53 ± 0.28 


27.69 ± 0.42 


21.93 ±0.35 


26.89 ± 0.53 


22.61 ± 0.35 


26.36 ± 0.44 


26.10 


4177 [Fe II] 


0.20 ± 0.05 


0.27 ± 0.07 












4233 Fe II 


0.18 ±0.04 


0.23 ± 0.05 


2.03 ± 0.23 


2.20 ± 0.25 








4244 [Fe II] 


0.29 ± 0.03 


0.38 ± 0.04 


4.32 ± 0.21 


4.66 ± 0.24 








4277 O II+[Fe II] 


0.19 ±0.02 


0.24 ± 0.03 


3.44 ± 0.21 


3.69 ± 0.23 








4287 [Fe II] 


0.48 ± 0.03 


0.61 ±0.04 


4.80 ± 0.20 


5.14 ±0.22 


0.24 ± 0.03 


0.26 ± 0.04 




4340 Hy 


39.69 ± 0.56 


49.92 ± 0.73 


40.96 ± 0.62 


46.08 ± 0.77 


42.59 ± 0.63 


47.00 ±0.71 


47.10 


4353 [Fe II] 


0.05 ± 0.02 


0.06 ± 0.03 


1.48 ±0.16 


1.56 ±0.17 








4359 [Fe II] 


0.53 ± 0.05 


0.66 ± 0.06 


4.24 ± 0.26 


4.49 ± 0.28 








4363 [O III] 


1.98 ±0.04 


2.45 ± 0.05 


2.80 ±0.31 


2.96 ± 0.34 


8.80 ±0.15 


9.59 ±0.17 




4387 He I 


0.30 ± 0.03 


0.36 ± 0.04 












4414 [Fe II]+0 II 


0.66 ± 0.06 


0.80 ± 0.07 


4.87 ± 0.52 


5.11 ±0.56 








4452 [Fe II] 


0.28 ± 0.05 


0.34 ± 0.06 


1.87 ±0.22 


1.95 ±0.23 








4458 [Fe II] 


0.20 ± 0.04 


0.23 ± 0.05 












4471 He I 


3.40 ± 0.06 


4.00 ± 0.07 


4.04 ± 0.63 


4.20 ± 0.68 


3.22 ± 0.09 


3.44 ±0.10 




4475 [Fe II] 


0.23 ± 0.05 


0.27 ± 0.06 


3.24 ± 0.30 


3.37 ± 0.32 








4493 [Fe II] 


0.20 ± 0.03 


0.23 ± 0.04 


1.06 ±0.15 


1.10 ±0.16 








4515 [Fen] 


0.13 ±0.02 


0.15 ±0.03 












4563 [Cr III]? 


0.09 ± 0.01 


0.10 ±0.02 












4571 Mg I 


0.15 ±0.02 


0.16 ±0.03 












4584 Fe II 


0.15 ±0.02 


0.17 ±0.02 


1.84 ±0.15 


1.88 ±0.16 








4597 [Co IV]? 


0.09 ± 0.02 


0.10 ±0.02 












4607 N II+[Fe III] 


0.14 ±0.02 


0.15 ±0.03 












4630 N II 


0.30 ± 0.03 


0.33 ± 0.04 


1.87 ±0.16 


1.90 ±0.16 








4640 N III 


0.19 ±0.04 


0.21 ±0.04 












4658 [Fe III] 


1.51 ±0.03 


1.64 ±0.04 


3.36 ±0.21 


3.39 ±0.21 


0.99 ± 0.05 


1.03 ±0.05 




4686 He II 


0.91 ±0.06 


0.98 ± 0.07 






2.15 ±0.06 


2.21 ±0.07 




4702 [Fe III] 


0.34 ± 0.03 


0.36 ± 0.04 






0.26 ± 0.04 


0.26 ± 0.04 




4712 [Ar IV] + He I 


0.67 ± 0.03 


0.71 ±0.04 






0.90 ± 0.07 


0.92 ± 0.07 




4728 [Fe II] 


0.16 ±0.03 


0.17 ±0.03 












4740 [Ar IV] 


0.17 ±0.01 


0.17 ±0.02 






0.83 ± 0.07 


0.84 ± 0.07 




4755 [Fe III] 


0.36 ± 0.03 


0.37 ± 0.03 












4770 [Fe III] 


0.14 ±0.02 


0.14 ±0.02 












4778 [Fe III] 


0.15 ±0.03 


0.15 ±0.03 












4814 [Fe II] 


0.27 ± 0.02 


0.28 ± 0.02 


2.06 ±0.16 


2.03 ±0.16 








4861 H/? 


100.00 ± 1.42 


100.00 ± 1.44 


100.00 ± 1.45 


100.00 ± 1.50 


100.00 ± 1.45 


100.00 ± 1.46 


100.00 


4881 [Fein] 


0.47 ± 0.02 


0.47 ± 0.02 






0.53 ± 0.04 


0.53 ± 0.04 




4890 [Fe II] 


0.16 ±0.01 


0.16 ±0.02 


1.14 ±0.10 


1.11 ±0.10 








4905 [Fe II] 


0.13 ±0.01 


0.13 ±0.01 


1.30 ±0.14 


1.26 ±0.14 


0.25 ± 0.02 


0.25 ± 0.02 
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4922 He I 


1.15 ±0.03 


1.12 ±0.03 






4930 [Fe III] 


0.20 ± 0.03 


0.19 ±0.03 






4943 O II 


0.09 ± 0.02 


0.09 ± 0.02 






4959 [0 III] 


124.43 ± 1.76 


119.38 ± 1.72 


102.66 ± 1.49 


99.02 ± 1.48 


4986 [Fe III] 


0.61 ± 0.02 


0.58 ± 0.02 


1.36 ±0.20 


1.31 ±0.20 


5007 [O III] 


394.68 ± 5.59 


371.56 ±5.34 


331.11 ±4.76 


317.19 ±4.68 


5016 He I 


1.38 ±0.03 


1.29 ±0.03 


2.32 ± 0.22 


2.22 ± 0.22 


5041 Si II 


0.25 ± 0.03 


0.23 ± 0.03 






5048 He I 


0.13 ±0.02 


0.12 ±0.02 






5056 Si II 


0.41 ± 0.03 


0.37 ± 0.03 






5159 [Fell] 


0.65 ± 0.02 


0.58 ± 0.02 


5.66 ±0.19 


5.31 ±0.18 


5169 Fell 


0.14 ±0.02 


0.12 ±0.01 


1.53 ±0.14 


1.43 ±0.13 


5199 [N I] 


1.79 ±0.03 


1.57 ±0.03 


4.83 ±0.17 


4.50 ±0.17 


5220 [Fe II] 


0.14 ±0.04 


0.12 ±0.03 






5235 Fe II? 


0.19 ±0.03 


0.17 ±0.02 


1.28 ±0.13 


1.19 ±0.13 


5262 [Fe II] 


0.29 ± 0.02 


0.25 ± 0.01 


1.68 ±0.12 


1.55 ±0.11 


5270 [Fe III] 


1.06 ±0.03 


0.90 ± 0.02 


2.82 ±0.12 


2.60 ±0.12 


5273 [Fe II] 


0.22 ± 0.01 


0.18 ±0.01 


2.52 ±0.11 


2.33 ±0.11 


5317 Fell? 


0.26 ± 0.02 


0.21 ± 0.01 


3.26 ±0.18 


2.99 ±0.17 


5334 [Fe II] 


0.21 ± 0.03 


0.17 ±0.03 


1.78 ±0.14 


1.63 ±0.13 


5363 [Ni IV] 


0.17 ±0.06 


0.14 ±0.05 


1.23 ±0.26 


1.12 ±0.24 


5376 [Fe II] 


0.20 ± 0.03 


0.17 ±0.02 


2.09 ± 0.29 


1.91 ±0.27 


5518 [CI III] 


0.38 ± 0.03 


0.29 ± 0.02 


0.86 ±0.13 


0.77 ±0.12 


5535Fell?Nll?Cll? 


0.12 ±0.08 


0.09 ± 0.06 






5538 [CI III] 


0.26 ± 0.02 


0.20 ± 0.01 






5755 [N n] 


1.62 ±0.03 


1.16 ±0.02 


6.04 ±0.18 


5.25 ±0.16 


5876 He I 


19.52 ±0.29 


13.50 ±0.21 


16.12 ±0.47 


13.81 ±0.42 


5958 Si II 


0.16 ±0.02 


0.08 ± 0.01 






5979 Si II 


0.36 ± 0.03 


0.24 ± 0.02 






6046 O I 


0.13 ±0.04 


0.09 ± 0.02 






6300 [O I] 


6.67 ±0.10 


4.06 ± 0.07 


9.55 ± 0.27 


7.81 ± 0.23 


6312 [S III] 


1.82 ±0.05 


1.10 ±0.03 


1.28 ±0.21 


1.05 ±0.18 


6347 Si II 


0.47 ± 0.02 


0.28 ± 0.02 






6364 [O I] 


2.36 ± 0.04 


1.41 ±0.02 


2.69 ±0.12 


2.18 ±0.10 


6365 [Ni II] 


0.16 ±0.02 


0.09 ± 0.01 






6371 Si II 


0.30 ± 0.02 


0.18 ±0.01 






6384 Fe II 


0.19 ±0.04 


0.11 ±0.02 






6456 Fe II 


0.25 ± 0.04 


0.15 ±0.02 


1.00 ±0.13 


0.81 ±0.10 


6516 Fell 


0.12 ±0.02 


0.07 ± 0.01 






6548 [N II] 


28.91 ± 0.41 


16.38 ± 0.26 


22.86 ± 0.38 


18.22 ± 0.33 


6563 Ha d 


514.82 ± 7.29 


290.71 ± 4.53 


358.69 d ±5.17 


286.36 d ± 4.60 


6583 [N II] 


86.96 ± 1.26 


48.81 ± 0.78 


72.57 ± 1.06 


57.64 ± 0.94 


6667 [Ni II] 


0.18 ±0.03 


0.10 ±0.02 






6678 He I 


6.36 ± 0.09 


3.48 ± 0.06 


3.50 ±0.10 


2.75 ± 0.08 


6716 [S II] 


28.86 ± 0.41 


15.63 ± 0.25 


37.16 ±0.55 


29.14 ±0.48 


6731 [S II] 


18.99 ±0.27 


10.25 ±0.16 


17.08 ± 0.27 


13.37 ± 0.24 


7002 O I 


0.09 ± 0.01 


0.04 ± 0.00 






7065 He I 


6.73 ±0.10 


3.34 ± 0.06 


3.98 ± 0.09 


3.02 ± 0.07 


7136 [Arm] 


14.37 ± 0.21 


7.01 ±0.12 


10.49 ± 0.26 


7.92 ±0.21 


7155 [Fe II] 


0.55 ± 0.03 


0.27 ± 0.01 


1.76 ±0.22 


1.32 ±0.17 


7172 [Fell] 


0.14 ±0.02 


0.07 ± 0.01 






7236 C II 


0.19 ±0.01 


0.09 ± 0.01 






7254 O I 


0.31 ±0.02 


0.15 ±0.01 






7281 He I 


0.94 ± 0.02 


0.44 ± 0.01 






7291 [Ca II] 


0.19 ±0.02 


0.09 ± 0.01 






7320 [0 II] 


5.32 ± 0.08 


2.49 ± 0.04 


3.47 ±0.10 


2.58 ± 0.08 


7330 [O II] 


4.38 ± 0.06 


2.04 ± 0.03 


2.35 ±0.11 


1.74 ±0.09 


7378 [Ni II] 


0.67 ± 0.02 


0.31 ±0.01 


1.66 ±0.12 


1.23 ±0.09 


7412 [Ni ill 


0.15 ±0.02 


0.07 ± 0.01 


0.44 ± 0.09 


0.32 ± 0.07 


7711 Fell 


0.08 ± 0.01 


0.03 ± 0.00 






7751 [Arm] 


3.56 ± 0.05 


1.52 ±0.03 


1.30 ±0.10 


0.94 ± 0.07 


7816 He I 


0.10 ±0.02 


0.04 ± 0.01 






8125 [Cm] 


0.23 ± 0.02 


0.09 ± 0.01 






8188 N I 


0.37 ± 0.03 


0.14 ±0.01 






8216 N I 


0.42 ± 0.06 


0.16 ±0.02 






8223 N I? 


0.41 ± 0.03 


0.16 ±0.01 






8229 Fe II 


0.26 ± 0.03 


0.10 ±0.01 






8242 N I 


0.27 ± 0.02 


0.10 ±0.01 






8299 P28 


0.09 ± 0.01 


0.03 ± 0.00 






8306 P27 


0.14 ±0.01 


0.05 ± 0.01 






8314 P26 


0.09 ± 0.01 


0.03 ± 0.00 






8323 P25 


0.16 ±0.01 


0.06 ± 0.00 






8334 P24 


0.31 ±0.02 


0.12 ±0.01 






8345 P23 


0.27 ± 0.02 


0.10 ±0.01 






8359 P22 


0.51 ±0.01 


0.19 ±0.01 






8374 P21 


0.37 ± 0.01 


0.14 ±0.00 






8392 P20 


0.40 ± 0.01 


0.15 ±0.00 






8413 P19 


0.45 ± 0.01 


0.17 ±0.00 







0.63 ± 0.05 



0.62 ± 0.05 



171.17 ±2.46 167.96 ±2.44 

1.25 ±0.05 1.22 ±0.05 

531.05 ±7.61 517.00 ±7.47 

1.78 ±0.06 1.73 ±0.06 



0.50 



0.58 ± 



13.16 ± 



0.05 



0.04 



0.47 



0.54 ± 



0.23 



11.27 ± 



0.05 



0.03 



0.20 



2.95 ± 0.09 2.39 ± 0.07 

1.62 ±0.08 1.32 ±0.06 



0.66 ±0.10 



0.53 ± 0.08 



2.91 ± 0.07 2.29 ± 0.06 

374.18 ±5.37 294.56 ± 4.62 

9.04 ±0.15 7.10 ±0.13 

3.90 ± 0.08 3.03 ± 0.07 

14.07 ±0.22 10.88 ±0.18 

10.99 ±0.17 8.49 ±0.15 



3.84 ± 0.07 
7.56 ±0.13 



0.68 ± 0.04 

2.42 ± 0.06 
2.08 ± 0.05 



2.87 ± 0.06 
5.60 ±0.11 



0.50 ± 0.03 

1.76 ±0.04 
1.51 ±0.04 



1.90 



0.05 



1.33 



0.04 



0.26 ± 0.03 
0.26 ± 0.04 
0.37 ± 0.05 
0.26 ± 0.04 
0.37 ± 0.04 
0.35 ± 0.04 



0.18 ±0.02 
0.18 ±0.03 
0.25 ± 0.03 
0.17 ±0.03 
0.25 ± 0.03 
0.24 ± 0.03 



282.00 



0.08 
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8438 P18 


0.64 ± 0.01 


0.24 ± 0.01 






0.37 ± 0.03 


0.25 ± 0.02 


0.31 


8446 O I 


3.93 ± 0.06 


1.47 ±0.03 


7.34 ±0.16 


5.03 ±0.12 


0.74 ± 0.04 


0.49 ± 0.03 




8470 P17 


0.62 ± 0.01 


0.23 ± 0.01 






0.36 ± 0.04 


0.24 ± 0.03 


0.36 


8490 Fe II 


0.19 ±0.01 


0.07 ± 0.01 












8505 P16 


0.76 ± 0.02 


0.28 ± 0.01 






0.41 ± 0.03 


0.27 ± 0.02 


0.44 


8548 P15 


0.90 ± 0.02 


0.33 ±0.01 






0.37 ± 0.04 


0.24 ± 0.02 


0.53 


8579? 


0.31 ±0.02 


0.11 ±0.01 












8601 P14 


1.25 ±0.02 


0.46 ± 0.01 






0.70 ± 0.04 


0.46 ± 0.03 


0.65 


8617 [Fe II] 


0.52 ± 0.02 


0.19 ±0.01 












8667 P13 


1.68 ±0.03 


0.60 ± 0.01 






0.76 ± 0.04 


0.50 ± 0.03 


0.81 


8683 N i 


0.48 ± 0.02 


0.17 ±0.01 












8753 P12 


2.19 ±0.04 


0.78 ± 0.02 


0.81 ±0.06 


0.55 ± 0.04 


1.23 ±0.08 


0.80 ± 0.05 


1.04 


8865 PI 1 


2.81 ±0.04 


0.98 ± 0.02 






1.62 ±0.05 


1.04 ±0.03 


1.35 


9017 P10 


3.66 ± 0.06 


1.24 ±0.02 






2.13 ±0.05 


1.35 ±0.04 


1.80 


9069 [S III] 


46.14 ±0.65 


15.58 ± 0.29 


21.14 ±0.33 


14.43 ± 0.29 


15.17 ±0.23 


9.68 ±0.19 




9095 Ca I] 


0.19 ±0.02 


0.06 ± 0.01 












9231 P9 


6.11 ±0.10 


2.01 ± 0.04 






1.42 e ±0.07 


0.89 e ± 0.04 


2.49 


9464 He I 


0.37 ± 0.02 


0.12 ±0.01 












9532 [S III] 


136.72 ± 1.94 


43.16 ±0.82 


59.08 ± 0.86 


38.54 ± 0.76 


41.96 ± 0.62 


25.96 ± 0.50 




9548 Pe 


10.08 ±0.15 


3.17 ±0.06 


3.26 ±0.12 


2.10 ±0.08 


3.51 ±0.09 


2.17 ±0.06 


3.57 


9850 [C I] 


0.57 ± 0.06 


0.17 ±0.02 



















b)NIR 


range 








10050 Pd 










9.56 ± 0.32 


5.75 ±0.21 


5.40 


10432 [Fe II]? 


0.18 ±0.06 


0.05 ± 0.02 












10502 Fe II 


0.40 ± 0.04 


0.11 ±0.01 


1.94 ±0.22 


1.20 ±0.14 








10551 ? 


0.27 ± 0.03 


0.07 ± 0.01 












10715 [Ni II] 


0.21 ±0.04 


0.06 ± 0.01 












10829 He I 


120.40 ± 1.71 


32.64 ± 0.66 


46.44 ± 0.76 


28.33 ± 0.63 


47.85 ± 2.23 


27.75 ± 1.36 




10910 He I 


1.58 ±0.11 


0.42 ± 0.03 












10941 Py 


20.45 e ± 0.32 


5.48 e ±0.12 


8.40 e ± 0.38 


5.10 c ±0.25 


13.90 ± 0.25 


8.03 ±0.18 


S.77 


11287 Ol+Fe II 


1.74 ±0.10 


0.45 ± 0.03 












11620 H 2 


0.18 ±0.04 


0.05 ± 0.01 












11850 H 2 


0.19 ±0.04 


0.05 ± 0.01 












11880 H 2 +[Fe II]? 


0.83 ± 0.05 


0.20 ± 0.01 












1 1970 He I 


0.51 ±0.05 


0.13 ±0.01 












12070 H 2 


0.09 ± 0.03 


0.02 ± 0.01 












12289 N I 


0.15 ±0.03 


0.04 ± 0.01 












12330 H 2 


0.52 ± 0.04 


0.13 ±0.01 












12380 H 2 +Fe II? 


0.19 ±0.03 


0.04 ± 0.01 












12525 ? 


0.90 ± 0.04 


0.21 ±0.01 












12570 [Fe II] 


6.14 ±0.11 


1.44 ±0.03 


4.77 ± 0.22 


2.77 ±0.14 








12790 He I 


3.83 ± 0.07 


0.89 ± 0.02 


2.77 ±0.19 


1.59 ±0.12 


0.88 ± 0.05 


0.48 ± 0.03 




12821 P/3 


62.09 ± 0.88 


14.36 ± 0.31 


30.54 ± 0.53 


17.59 ± 0.42 


27.22 ± 0.42 


14.78 ± 0.33 


15.70 


13165 O I 


0.92 ±0.10 


0.21 ±0.02 












14967 Br25 


0.34 ± 0.03 


0.07 ± 0.01 








... 


3.06 


15005 Br24 


0.28 ± 0.03 


0.06 ± 0.01 








... 


3.06 


15043 Br23 


0.25 ± 0.02 


0.05 ± 0.00 








... 


).07 


15087 Br22 


0.34 ± 0.03 


0.07 ± 0.01 








... 


).08 


15137 Br21 


0.35 ± 0.03 


0.07 ± 0.01 








... 


).09 


15196 Br20 


0.50 ± 0.03 


0.10 ±0.01 








... 


).ll 


15265 Brl9 


0.40 ± 0.02 


0.08 ± 0.01 








... 


3.12 


15346 Brl8 


1.19 ±0.04 


0.24 ± 0.01 








... 


).15 


15443 Brl7 


1.12 ±0.05 


0.23 ± 0.01 








... 


).18 


15530 ? 


0.33 ± 0.03 


0.07 ± 0.01 












15561 Brl6 


1.05 ±0.07 


0.21 ± 0.01 








... 


).21 


15705 Brl5 


1.29 ±0.05 


0.26 ± 0.01 








... 


3.26 


15885 Brl4 


1.49 ±0.04 


0.30 ± 0.01 






0.44 ± 0.05 


0.22 ± 0.03 


).31 


16010 H 2 


0.10 ±0.03 


0.02 ± 0.01 












16114 Brl3 


2.08 ± 0.04 


0.42 ±0.01 






0.71 ± 0.04 


0.36 ± 0.02 


).39 


16410 Brl2 


2.24 ± 0.05 


0.44 ±0.01 






0.77 ± 0.04 


0.39 ± 0.02 


3.50 


16436 [Fe II] 


6.28 ±0.10 


1.25 ±0.03 


3.14 ±0.18 


1.71 ±0.11 








16640 [Fe II] 


0.22 ± 0.05 


0.04 ± 0.01 












16769 [Fe II] 


0.36 ± 0.03 


0.07 ± 0.01 












16811 Brll 


3.51 ± 0.06 


0.69 ± 0.02 






1.02 ±0.08 


0.52 ± 0.04 


).65 


16873 Fe II 


0.40 ± 0.03 


0.08 ± 0.01 












17006 He I 


1.16 ±0.03 


0.23 ± 0.01 






0.55 ± 0.04 


0.28 ± 0.02 




17367 BrlO 


5.10 ±0.08 


0.99 ± 0.02 


2.60 ± 0.37 


1.41 ± 0.20 


1.92 ±0.06 


0.97 ± 0.03 


).87 


18179 Br9 


744 e ±0.31 


1.42 e ±0.06 










1.21 


18693 ? 


9.90 ± 0.34 


1.88 ±0.07 












18756 Pa 


42.33 e ± 0.77 


8.02 e ± 0.21 






13.11 e ± 0.23 


6.55 e ±0.16 


31.90 


19450 Br<S 


10.06 ±0.17 


1.89 ±0.05 


2.63 ± 0.22 


1.40 ±0.12 


1.18 ±0.09 


0.59 ± 0.05 


1.73 


19540 He l+[Fe II] 


0.73 ± 0.07 


0.14 ±0.01 












19570 H 2 


0.61 ± 0.03 


0.11 ±0.01 






0.58 ± 0.07 


0.29 ± 0.03 




20340 H 2 


0.40 ± 0.08 


0.07 ± 0.02 












20586 He I 


8.22 ±0.13 


1.52 ±0.04 


1.48 ±0.14 


0.79 ± 0.08 


2.02 ± 0.09 


1.00 ±0.05 




20730 H 2 


0.33 ± 0.03 


0.06 ± 0.01 
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20888 Fen 0.19 ±0.03 0.04 ± 0.01 

21130 He I 0.57 ±0.04 0.11 ±0.01 

21220 H 2 1.21 ±0.03 0.22 ± 0.01 ... ... 0.56 ±0.06 0.28 ±0.03 ... 

21620 He I 0.42 ± 0.04 0.08 ± 0.01 

21661 Bry 16.73 ± 0.24 3.07 ± 0.07 5.38 ±0.19 2.85 ±0.12 5.03 ±0.12 2.48 ± 0.07 2.63 

22181? 0.37 ±0.04 0.07 ± 0.01 

22230 H 2 0.26 ±0.05 0.05 ± 0.01 

22480 H 2 0.20 ±0.04 0.04 ± 0.01 

c) MIR range f 

10.51//m[S IV] 
n.SljumlNe II] 
15.55/jm[Ne III] 
18.71//m[S III] 

C(H/3) 0.74 0.27 0.31 

FiUfif 344.20 63.59 121.40 

EW(abs) O20 035 (L15 

a F(A) is observed flux, F(A)/F(HfS) rado designadons with X100. 
b 1(A) is extinction-corrected flux, /(/l)//(H/3) ratio d esignations with xlOO. 
c hydrogen recombinationdine relative intensities by Hummer & Storey 1 1987). 

d The observed flux of Ho- in Haro 11C is F(Ho-)=3.350xl0~ 14 erg s cirT 2 . In this Table we show the F(Ha)/F(Hfl) ratio corresponding to F(Ha')=2.281xlO~ 
erg s _1 cm -2 obtained after subtraction of the P Cyg component. 
c affected by the telluric absorption. 
f F(A)/F(Hfi) ratio from lWuet al.li2008b . 
1 erg s - ' cm -2 . 



12.76 


±0.10 


9.07 


±0.14 


28.51 


±0.40 


13.14 


±0.12 




0.74 




344.20 




0.20 



